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ABSTRACT 

We present the results of a combined investigation of the spectral and kHz QPO evolution around the Z-track in GX5-1 based 
on high-quality Rossi-XTE data. In spectral analysis, we find that the Extended ADC emission model provides very good fits 
to all of the spectra, and the results point clearly to a model for the nature of the Z-track in this source, in agreement with 
previous results for the similar source GX 340+0. In this model, at the soft apex of the Z-track, the mass accretion rate M 
is at its minimum and the neutron star has its lowest temperature; but as the source moves along the normal branch, the 
luminosity of the Comptonized emission increases, indicating that M increases and the neutron star gets hotter. The measured 
flux / of the neutron star emission increases by a factor of ten becoming super-Eddington, and we propose that this causes 
disruption of the inner disk and the formation of jets. In flaring, the luminosity of the dominant Comptonized emission from 
the accretion disk corona is constant, while the neutron star emission increases, and we propose for the first time that flaring 
consists of unstable nuclear burning on the neutron star, supported by the agreement between the measured mass accretion rate 
per unit area rh at the onset of flaring and the theoretical critical value at which burning becomes unstable. There is a striking 
correlation between the frequencies of the kHz QPO and the ratio of the flux to the Eddington value: ///Edd, suggesting an 
explanation of the higher frequency QPO and of its variation along the Z-track. It is well known that a Keplerian orbit in the 
disk at this frequency corresponds to a position some distance from the neutron star; we propose that the oscillation always 
occurs at the inner disk edge, which moves radially outwards on the upper normal and horizontal branches as the measured 
increasing radiation pressure increasingly disrupts the inner disk. 

Key words. Accretion: accretion disks - acceleration of particles - binaries: close - line: formation - stars: neutron - X-rays: 
binaries - X-rays: individual (GX5-1) 



1. Introduction 

The Z-track sources form the brightest group of low mass 
X-ray binaries (LMXB) containing a neutron star, with lu- 
minosities at or above the Eddington limit. They are char- 
acterised by having three distinct branches in a hardness- 
intensity diagram: the horizontal branch (HB), the nor- 
mal branch (NB) and the flaring branch (FB) (Hasinger 
& van der Klis 1989) showing that major physical changes 
take place within the sources but the nature of these has 
not been understood. It has been widely thought that 
the physical and spectral changes are driven by change 
of a single parameter along the Z-track, presumably the 
mass accretion rate (Priedhorsky et al. 1986), assumed 
to increase monotonically in the direction HB - NB - FB. 
However, as discussed below, the evidence for this is rather 
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limited. Moreover the variation of the X-ray intensity does 
not obviously support this as the intensity does not in- 
crease monotonically in the direction HB - NB - FB, but 
decreases on the normal branch. Arguably the most impor- 
tant feature of the Z-track sources is the detection of radio 
emission showing jets to be present but only in the upper 
normal and horizontal branches (Penninx 1989). The pres- 
ence of jets was dramatically demonstrated by extended 
radio observations of the Z-track source Sco X-l (Fomalont 
et al. 2001) which revealed radio condensations moving 
away from the source with velocity v/c of 0.45. Thus the Z- 
track sources uniquely provide the possibility of determin- 
ing the physical conditions within the sources on the part 
of the Z-track where radio is detected so telling us the con- 
ditions needed for the launching of jets. Apart from this 
an understanding of the Z-track sources is essential to the 
basic understanding of LMXB in general. Extensive work 
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has been carried out on the timing properties of Z-track 
sources (van der Klis et al. 1987; Hasinger & van der Klis, 
1989) and revealed the existence of quasi-periodic oscil- 
lations (QPO) which change along the Z-track. However, 
analysis of the timing properties has not provided an ex- 
planation of the Z-track phenomenon. Spectral analysis is 
more likely to reveal the nature of the physical changes 
taking place as directly showing changes in the emission 
components during the spectral evolution along the Z- 
track, but spectral studies of the sources have been hin- 
dered by lack of agreement over the emission model to be 
used. 

1.1. Spectral studies of the Z-track sources 

The nature of X-ray emission in LMXB has been con- 
troversial essentially because two radically different types 
of model for the continuum are capable of describing the 
spectra in which a dominant Comptonized component and 
also a thermal component are clearly present. Analysis of 
the dipping class of LMXB led to the "Extended ADC" 
model in which the dominant component is Comptonized 
emission from the accretion disk corona (ADC) plus ther- 
mal emission from the neutron star (Church & Baluciflska- 
Church 1995, 2004). Both spectral and timing analysis 
showed the Comptonizing ADC to be a region extend- 
ing over a substantial fraction of the inner accretion disk 
(typically 15%) and thin (H/r « 1) as discussed in more 
detail below. The other model is the Eastern model com- 
prising multi-colour blackbody emission from the inner 
accretion disk plus Comptonized emission from a central 
region close to the neutron star (Mitsuda et al. 1989). 
Previous spectral fitting of the Z-track sources has mostly 
been based on use of the Eastern model as summarized 
below. However, there is now a body of evidence that the 
Comptonizing region is extended (Sect. 1.2) and is not a 
small central region. Moreover, when the Eastern model 
was applied to the Z-track sources, it was difficult to in- 
terpret the variation of spectral parameters and a physical 
explanation was often not given. 

An extensive investigation of spectral evolution along 
the Z-track in several sources using Exosat data was car- 
ried out by Schulz et al. (1989) and by Schulz and Wijers 
(1993) using a model with neutron star blackbody emis- 
sion plus Comptonization. Hasinger et al. (1990) fitted 
Ginga data on Cyg X-2 with both Western and Eastern 
models. In the Western model the main change was an 
increase of blackbody temperature on climbing the nor- 
mal branch. In the Eastern model, there was a systematic 
change in the Comptonization y-parameter on the nor- 
mal branch. Because of the ability of both models to fit, 
it was not thought possible to distinguish between phys- 
ical models. Consideration of multi-wavelength data and 
of the timing properties of the source led the authors to 
conclude that the mass accretion rate M increased mono- 
tonically in the direction HB - NB - FB although X-ray 
intensity decreased on the NB. A possible explanation was 



that a thickening of the inner disk obscured the X-ray 
emission. The conclusion about M led Hasinger et al., fol- 
lowing Hasinger (1988) and Lamb (1989), to propose a 
scenario in which M was small on the HB so that a thin 
disk existed together with a magnetosphcrc in which ra- 
dio emission was generated. Increasing M resulted in the 
magnetosphere being enveloped by a vertically expanding 
disk at the hard apex, suppressing radio emission. At the 
soft apex M exceeds the Eddington value so that strong 
radiation pressure led to mass ejection above and below 
the disk. Asai et al. (1994) fitted Ginga data on the NB 
and FB spectra of GX5-1 with both Western and Eastern 
models, but using the Eastern model in its original form 
consisting of disk blackbody plus blackbody from the neu- 
tron star. This work concentrated on the emission feature 
found at 10 keV. 

Spectral fitting of Cyg X-2 was carried out by Done 
et al. (2002) using the Eastern model assuming that the 
seed photons for Comptonization come from the neu- 
tron star or inner disk and so should be described by a 
blackbody spectrum. It was also assumed that the central 
Comptonizing cloud would illuminate the disk so that re- 
flection should be present (although in the Z-track sources 
the height of the inner disk is much greater than that 
of the neutron star). The spectral fitting indicated an 
increasing flux of the disk blackbody in the lower nor- 
mal and flaring branches where it was assumed that M 
increased; however, the disk temperature appeared con- 
stant so it was argued that the inner disk radius must 
increase. The Comptonized flux also remained constant 
when it might be expected to increase and it was not 
clear what caused the spectral parameters to vary in the 
way found. Agrawal & Sreekumar (2003) used the Eastern 
model to fit Rossi-XTE spectra of GX 349+2 in the NB 
and FB using the form diskbb + COMPTT, i.e. disk black- 
body plus Comptonization with the seed photons having 
a Wien distribution. The inner disk temperature first in- 
creased on the NB then decreased in flaring; the inner disk 
radius first decreased then increased, suggesting that the 
inner disk edge was moving to larger radial positions, and 
the authors proposed that there was thus an outflow of 
material . The electron temperature of the Comptonizing 
plasma first increased on the NB then decreased in flaring 
while the optical depth of the region first decreased then 
increased. It was suggested that the initial heating caused 
a density decrease and reduction in optical depth; on the 
FB the outflow proposed added density to the hot cen- 
tral corona causing the increase of optical depth on this 
branch. 

Spectral fitting of broadband BeppoSAX data was car- 
ried out for GX 17+2, GX 349+2 and Cyg X-2 (di Salvo et 
al. 2000, 20001, 2002). In the case of CygX-2, the Eastern 
model was applied together with a power law to repre- 
sent a high energy tail, and it was argued that all three 
sources could be fitted by this model. In Cyg X-2 spectral 
evolution consisted of an increase of disk blackbody tem- 
perature and decrease of blackbody radius moving from 
the HB to the NB, suggesting that the inner disk ra- 
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clius was shrinking although an explanation of this was 
not proposed. There were changes in temperature of the 
Comptonizing plasma and on the normal branch, the in- 
ferred radius of this region fell to 2 - 3 km, incompati- 
ble with the neutron star being the source of seed pho- 
tons unless there was a marked change in geometry from 
the normally assumed spherical Comptonizing cloud of the 
Eastern model. D'Ai et al. (2007) investigated the Z-track 
in ScoX-1 using extensive numbers of observations with 
RXTE. The Eastern model in the form used by di Salvo et 
al. was used. The disk blackbody temperature was highest 
on the flaring branch as expected on the standard assump- 
tion that M is largest here. The electron temperature of 
the Comptonizing plasma increased around the Z having 
its highest values in flaring, while at the same time the 
optical depth decreased; however, the reasons for these 
changes were not clear. 

It is difficult to find consistency in the spectral fit- 
ting results and with the difficulty of interpreting the 
evolution of spectral parameters along the Z-track it is 
fair to say that there is no general concensus on the na- 
ture of the Z-track. Much effort was devoted to producing 
a theoretical model of of accretion in the inner disk in 
LMXB which proposed detailed physical changes taking 
place along the Z-track (Lamb 1989; Miller 1990; Psaltis 
et al. 1995). This involves a magnetosphcre at the inner 
accretion disk, and changes in the geometry and extent of 
the magnetosphcre. The model assumes the main element 
of the Eastern model that is, a small central Comptonizing 
region. In the present work we take the approach that the 
evidence favours an extended Comptonizing region and 
we test the hypothesis that the Extended ADC model may 
provide a physically reasonable and consistent explanation 
of the Z-track source GX5-1. 

1.2. The extended ADC 

Our work on the dipping class of LMXB previously 
led to our proposal of the extended ADC model com- 
prising blackbody emission from the neutron star plus 
Comptonizcd emission from an extended, thin ADC ex- 
isting as a hot layer above the accretion disk. In the dip- 
ping sources with inclination angles between 65° and 85°, 
reductions in X-ray intensity take place at the orbital pe- 
riod generally thought to involve absorption in the bulge in 
the outer disk where the accretion flow impacts. Spectral 
evolution cannot be described in terms of absorption of a 
single emission component, but is more complex strongly 
constraining allowable models. It provides direct evidence 
for the extended nature of the Comptonizcd emission as 
this spectral component is removed only slowly during dip- 
ping showing the emitter to be extended. Secondly, the 
ADC size can be measured directly by the technique of dip 
ingress timing since the ingress time depends on the size 
of the major emission component: the Comptonized emis- 
sion. This has shown that the radial extent of the ADC is 
typically 50 000 km (Church & Balucinska-Church 2004) 



also showing that the ADC, like the disk is in general 
thin, with height H at radial position r having H/r << 
1. The extended ADC model was shown to give very good 
fits to the complex spectral evolution in dipping in many 
observations of the dipping LMXB (Church et al. 1997, 
1998a,b, 2005; Balucihska-Church et al. 1999, 2000; Smale 
et al. 2001; Barnard et al. 2001). In addition, it was shown 
able to fit the spectra of all types of LMXB in a survey of 
LMXB made with ASCA (Church & Balucihska-Church 
2001) and so was able to fit sources of all inclination an- 
gles. 

The evidence favouring the extended ADC is thus 
strong and it is difficult to avoid the conclusion that the 
main assumption of the Eastern model is invalid. 

The main difference between the spectral forms of the 
extended ADC model and the Eastern model relates to 
the seed photons. In the Eastern model it is usually as- 
sumed that these originate on the neutron star and can 
be described by a single temperature blackbody or by the 
Wien approximation. In the extended ADC model the ma- 
jor source of seed photons must be the thermal emission of 
the disk below the ADC. Because of the temperature gra- 
dient T{r) in the disk, the outer portions of the ADC will 
be above relatively cool parts of the disk providing a large 
population of soft seed photons (kT ~ 0.1 keV or less) and 
the seed photon spectrum is very different from a black- 
body at ^1 keV. Comptonization is well represented by a 
power law extending to energies lower than 0.1 keV, with 
a cut-off at high energies related to the maximum electron 
energy (see Church & Balucinska-Church 2004 for a full 
discussion). The use of a simple blackbody for the neu- 
tron star emission in applying the extended ADC model 
reflects the general lack of evidence for modification of the 
blackbody in the atmosphere of the neutron star (Church 
et al. 2002). Taking into account many spectral studies 
of LMXB, including the spectra of X-ray bursts, the only 
evidence for modification is in a minority of bursts, with 
no evidence in non-burst emission. Church et al. showed 
that whether modification occurs depends critically on the 
electron density in the neutron star atmosphere and that 
this is poorly known. 

1.3. A model for the Z-track 

We previously applied the extended ADC model to high 
quality Rossi- X-ray Timing Explorer data from the Z- 
track source GX 340+0 (Church et al. 2006), and found 
that it fitted the spectra well at all positions on the Z- 
track, and also clearly suggested a physical model for the 
phenomenon. In this model, the soft apex between NB and 
FB is a state of the source in which the intensity and the 
mass accretion rate were low. Moving onto the NB, there 
was a large increase in X-ray intensity (count s _1 ) and an 
increase in the total broadband luminosity by about a fac- 
tor of two caused by an increase in the luminosity Ladc 
of the Comptonized emission of the ADC. It was argued 
that this increase indicated that M was increasing on the 
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NB and was supported by a marked increase in the black- 
body temperature of the neutron star from ~1 keV to ^2 
keV and even higher values on the HB. This heating of 
the neutron star means a large increase in T 4 and so the 
radiation pressure became very high. There was a simul- 
taneous decrease in emitting area on the neutron star. In 
this situation, it is appropriate to compare the flux emit- 
ted per unit area of the star with the Eddington value 
(see Sect. 4.1), and the flux increased from a low value 
till on the HB, it was three times super-Eddington. It was 
argued that this would disrupt the inner disk deflecting 
the accretion flow into the vertical direction leading to jet 
formation, and explaining why jets are detected on the 
upper NB and HB where the radiation pressure was high. 
In flaring it was found that Labc was constant implying 
that M was constant although the neutron star luminos- 
ity increased. The mass accretion rate per unit area of the 
neutron star closely agreed with the theoretical value for 
the onset of unstable nuclear burning, and it was proposed 
thus flaring was due to this. 

In the present work, we apply the extended ADC 
model to high-quality Rossi-XTE observations of the Z- 
track source GX5-1 to test the hypothesis that the same 
physical explanation can explain a second Z-track source. 
Spectral and timing studies were made using the same 
selections of data allowing the evolution of kHz QPO fre- 
quency to be directly correlated with the spectral evo- 
lution. This revealed a striking dependence of kHz QPO 
frequency on the neutron star radiation pressure. In the 
next section kHz QPO are briefly reviewed. 

1.4. The kHz QPO 

Quasi periodic oscillations (QPO) have been the subject 
of intense study for many years, and a major discov- 
ery of the Rossi X-ray Timing Explorer was that of the 
twin kilohertz QPO in LMXB. Predictions of millisec- 
ond variability pre-date this considerably and Sunyaev 
(1973) argued that hot clumps of material orbiting in 
the inner disk around a black hole or neutron star would 
lead to quasiperiodic variability on millisecond timescales. 
Kilohertz QPO have now been detected in many LMXB 
(see reviews of van der Klis 2000, 2006), most of these dis- 
playing two peaks at frequencies such as at 600 and 900 
Hz, and in the Z-track sources both frequencies increase 
during movement along the HB towards the hard apex. 
Rather similar QPO have also been seen in a number of 
black hole binaries. 

In almost all QPO models the higher frequency peak 
is associated with the orbital frequency at some distance 
into the accretion disk away from the neutron star; for 
example, a frequency of 900 Hz corresponds to a radial 
position r of about 18 km. In the sonic point beat fre- 
quency model (Miller et al. 1998) radiation drag reduces 
the angular momentum in the disk within a few stellar 
radii of the neutron star leading to a sonic point where 
the radial velocity becomes high (Lamb 1989; Miller 1991; 



Miller & Lamb 1993, 1996). The upper frequency v 2 is as- 
sociated with inhomogencitics orbiting at the sonic point 
in the inner disk, which are irradiated by X-rays from hot 
regions rotating on the neutron star leading to a beat- 
ing of the upper frequency with the spin frequency ^ S pin 
so producing the lower frequency QPO. Although earlier 
observations found the frequency difference Av constant 
within errors consistent with the model, improved data 
showed that Av varied around the Z-track. The detection 
of twin kHz QPO and X-ray burst oscillations in the ac- 
creting millisecond pulsar SAX J1808. 4-3658 (Wijnands et 
al. 2003) led to the conclusion that Av ~ v sp i n in some 
sources, but ^ sp in/2 in other sources, inconsistent with the 
model, leading Lamb & Miller (2003) to propose a rela- 
tivistic disk-spin resonance model. 

The relativistic precession model (Stella & Vietri 1998) 
identifies v 2 with motion of inhomogencitics at the inner 
disk edge and the lower frequency v\ with precession of the 
inhomogeneities. The disk resonance model (Abramowicz 
& Kluzniak 2001) recognizes the preference of the kHz 
QPO to have a ratio of 3:2 in their frequencies suggest- 
ing a non-linear resonance between normal modes of the 
disk. This model predicted that twin peaks should also 
exist in black hole binaries, and the discovery of a sec- 
ond peak in GRO J1655-40 (Strohmayer 2001) was a suc- 
cess of the model. Comparison with data (Abramowicz et 
al. 2003; Kluzniak et al. 2007) for both neutron star and 
black hole systems supports the model: black hole sys- 
tems accurately follow the 3:2 ratio; LMXB show some 
scatter but a histogram of the ratio peaks strongly at 3:2. 
However, Belloni et al. (2005) have argued that in the case 
of ScoX-1, the peak in the distribution at a ratio of 1.5 
was only marginally significant. 

1.5. The assumed monotonic variation of mass 
accretion rate 

The assumption of M increasing in the direction HB - NB 
- FB generally, but not universally made, is not supported 
by our previous results for GX 340+0, and so we briefly 
review evidence for and against this. 

The evidence for this was based originally on a multi- 
wavelength campaign on CygX-2 (Hasinger et al. 1990) in- 
volving X-rays, ultraviolet, optical and radio data. Results 
from IUE (Vrtilek et al. 1990) had increasing intensity for 
Z-track movement HB - NB - FB implying an increase of 
M\ however, strong variability interpreted as flaring may 
have been due to X-ray dipping in which case there is no 
strong correlation of UV with track position. 

The properties of most power spectral features corre- 
late well with position on the Z-track; for example, the 
frequencies of the HBO (horizontal branch oscillations), 
i.e. low frequency QPO, increase around the Z-track in 
the direction HB - NB from ~5 to ^60 Hz (e.g. Homan et 
al. 2002). Similarly, NBO (normal branch oscillations) are 
seen on the NB and in two sources, ScoX-1 and GX 17+2, 
oscillations are seen in the FB, but only close to the soft 
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apex, and there is a continuous increase of QPO frequency 
from lower NB to FB (Casella et al. 2006), i.e. in crossing 
the soft apex between these branches. Similarly, the fre- 
quencies of kilohertz QPO increase along the HB where 
they are observed (e.g. van der Klis 2000). It is general 
practice to display QPO properties as a function of arc 
distance S measured along the Z-track in the direction 
HB - NB - FB, and this shows that QPO properties are 
generally well-correlated with S (e.g. Dieters & van der 
Klis 2000) implying a dependence on a single parameter, 
i.e. M. It is also generally though that the soft apex repre- 
sents the position at which M increases to the Eddington 
value (Lamb 1989). In addition, the frequencies of HBO 
and kHz QPO in many models are directly related to M 
(see Homan et al. 2002). 

However, in recent years, some authors are beginning 
to doubt the general assumption of a monotonic increase 
in M. Some results conflict with the assumption: e.g. in 
GX 17+2, the frequency of the HBO first increases, but 
then decreases on the normal branch (Wijnands et al. 
1996; Homan et al. 2002). In ScoX-1 in the neighbour- 
hood of the soft apex there is a rapid increase of QPO fre- 
quency (Casella et al. 2000; Dieters & van der Klis 2000) 
with S which is inconsistent with the assumption that £ is 
a measure of M as the change is over a very short part of 
the Z-track. Moreover the approximate constancy of NBO 
frequency in general over most of the NB does not indicate 
a dependency on M. Also, in GX17+2, the properties of 
the X-ray bursts sometimes observed do not depend on 
Z-track position (Kuulkers et al. 2002). 

In our previous work on GX 340+0 (Church et al. 2006) 
we proposed that M increased between the soft apex and 
the hard apex, while on the FB, M was constant. The X- 
ray intensity, the measured broadband luminosity of the 
source and the luminosity of the Comptonized emission 
all increased substantially on the NB strongly suggesting 
that M is not decreasing. The argument is often made 
that the X-ray intensity cannot be a good indicator of S 
and thus of M (e.g. Dieters & van der Klis 2000) because 
the intensity increases on both the NB and FB, moving 
away from the soft apex. However, our work on GX 340+0 
shows this not to be a secure argument, since the intensity 
increases on the NB because M increases, while on the FB, 
M is constant, but strong nuclear burning on the neutron 
star causes the intensity increase. One aim of the present 
work is to test whether similar results are found in GX 5- 1 . 

1.6. GX5-1 

GX5-1 is the second brightest persistent LMXB (after 
ScoX-1), discovered in 1968 (Fisher et al. 1968; Bradt 
et al. 1968) as a bright X-ray source in the Sagittarius 
region and classified as a Z-track source by Hasinger & 
van der Klis (1989). A three-year study of the source was 
made using Ginga (van der Klis et al. 1991) which re- 
vealed an extended horizontal branch and an upcurving 
at the end suggesting a change of behaviour, effectively 



a fourth branch (Lewin et al. 1992). Analysis of Exosat 
and Ginga data revealed flaring branch data for the first 
time (Kuulkers et al. 1994) and a detailed investigation 
of correlated spectral and timing behaviour was made. 
It was the first Z-track source found to have QPO (van 
der Klis et al. 1985) and kilohertz QPO were detected 
by Wijnands et al. (1998). Radio emission has been de- 
tected from GX5-1, e.g. by Penninx (1989) implying the 
existence of jets, and it is known to be strongest on the 
horizontal branch (Bcrendsen et al. 2000), as for other Z- 
track sources. Analysis of 1987 Ginga data (Asai et al. 
1994) and ASCA data (Asai et al. 2000) showed there was 
no evidence for an iron K line, but a Gaussian feature ex- 
isted at 10 keV. To date no satisfactory explanation exists 
for this line. 

2. Observations and analysis 

GX5-1 has been observed several times with RXTE. We 
examined the data in the heasarc archive remotely pro- 
ducing hardness-intensity plots and selected data which 
spanned the full Z-track in a relatively short time. It is 
important not to use data in which sideways movement 
of the Z-track takes place, broadening the track or even 
giving several shifted Z-tracks, as we are interested in the 
changes taking place around a single Z-track. Similarly 
data separated by months should not be co-added on the 
assumption that the source was in the same state because 
portions of the Z-track approximately overlap. Thus a sin- 
gle observation including all branches is preferred, and in 
view of the very high count rates in the RXTE instru- 
ments, there is no need to add further data as the count 
statistics are already excellent. We thus selected the obser- 
vation of 1998 November 21 - 22 having the best coverage 
of the Z-track in a single observation, consisting of four 
sub-observations spanning 90 ksec. 

Data from both the Proportional Counter Array 
(PCA, Jahoda et al. 1996) and the high-energy X-ray tim- 
ing experiment (HEXTE) were used. The PCA was in 
Standard2 mode with 16s resolution and examination of 
the housekeeping data revealed that all five of the xenon 
proportional counter units were on during the observa- 
tion. Standard screening criteria based on the housekeep- 
ing data were applied to select data having an offset be- 
tween the source and telescope pointing of less than 0.02° 
and elevation above the Earth's limb greater than 10°. 
Data were extracted from the top layer of the detector us- 
ing both left and right anodes. Analysis was carried out us- 
ing the FTOOLS 5.3.1 package. Lightcurves were generated 
in the band 1.9 - 18.5 keV and background files for each 
PCA data file produced using the facility pcabackest, 
applying the latest "bright" background model for Epoch 
3 of the mission (1996 April 15 - 1999 March 22) which 
includes the present observation. Deadtime correction was 
carried out on both source and background files prior to 
subtraction. Light curves were also made in sub-bands of 
the above range, and a hardness ratio defined as the ratio 
of the intensity in the band 7.3 - 18.5 keV to that in the 
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Fig. 1. Background-subtracted deadtime-corrected PCA lightcurve of the November 1998 observation of GX5-1 with 64 s 
binning. 



4.1 - 7.3 keV band. Tests were made on suitable binning 
for the lightcurves and Z-track. In order to reduce scat- 
ter in the plot of hardness versus intensity and produce a 
relatively narrow, well-defined Z-track, it was found that 
a binning of 176 s was optimum. In Fig. 1 we show the 
background-subtracted, deadtime-corrected lightcurve in 
the band 1.9 - 18.5 keV with the binning reduced from 
176 to 64 s to allow more detail to be seen. In Fig. 2 we 
show the background-subtracted, deadtime corrected Z- 
track with 176 s binning. We identified each subsection of 
the lightcurve with its position on the Z-track by making 
a hardness-intensity plot for each and finding where this 
lay on the overall Z-track. Each of the 14 subsections was 
selected using the ftool maketime. It can be seen that 
the intensity of the source varied between two well-defined 
limits corresponding to the Soft Apex and the Hard Apex 
of the Z-track as discussed below. 

HEXTE data were also extracted as lightcurves and 
spectra from Cluster 1 of the HEXTE instrument using 
the ftool hxtlcurve which also provided background 
files. Deadtime correction was applied using the deadtime 
coefficients file of February 2000. This allowed simultane- 
ous fitting of PCA and HEXTE spectra, the wider energy 
band resulting in spectral fitting being better constrained 
that with the PCA alone. 
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Fig. 2. Z-track of the observation derived from background- 
subtracted and deadtime-corrected lightcurves with 176 s bin- 
ning. The boxes show the ranges of intensity and hardness ratio 
used for the selection of PCA and HEXTE spectra, and timing 
data. 



3. Results 

From Fig. 1 the strong variability of GX5-1 during these 
observations can be seen, and the lightcurve is labelled 
to show the identifications with the branches of the Z- 
track. At the start of the observation there is strong flar- 
ing with individual flares lasting several thousand seconds. 
The source then moves to the normal branch with the in- 
tensity increasing. Following movement along the normal 
branch to the hard apex, the source reaches the horizon- 
tal branch, the intensity then falling on this branch to 
intermediate values. 



Figure 2 shows the corresponding Z-track having full 
coverage of the Z except for the lower part of the normal 
branch, as the source was in this position during a data 
gap between sub-observations. However, this gap was not 
a problem as it can be seen from the results in Figs. 4, 5 
and 6 that there is no ambiguity in the behaviour of the 
source in this part of the Z-track. In the present work, 
we carry out both spectral fitting and timing analysis for 
data selected at a sequence of positions along the Z-track. 
This approach thus allows direct comparison of the QPO 
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Table 1. Counts accumulated in the PC A and HEXTE for 
each selection; the selections are numbered 1 to 11 starting 
from the HB end of the Z-track 



selection 


PCA Count 


HEXTE Count 


1 


2.04xl0 6 


2. 98 xlO 3 


2 


4.30 xlO 6 


6.38xl0 3 


3 


10.23xl0 6 


11.69xl0 3 


4 


4.69 xlO 6 


6.66xl0 3 


5 


15.71xl0 6 


15.16 xlO 3 


6 


2.17xl0 6 


3.81xl0 3 


7 


9.78 x10 s 


10.30xl0 3 


8 


6.78 x10 s 


6.17xl0 3 


9 


3.05 x10 s 


4.52xl0 3 


10 


1.63 x10 s 


2.08xl0 3 


11 


1.74 x10 s 


1.89xl0 3 



results with spectral fitting results for exactly the same 
data. 

Our previous work has suggested the importance of 
selecting data lying on a smooth curve along the cen- 
tre of the Z-track (Barnard et al. 2003) so that spectral 
changes reflect changes along the track and not perpen- 
dicular to it. PCA spectra were generated about equally- 
spaced along the Z-track selecting data within boxes in 
hardness-intensity by use of a good time interval (GTI) 
file for each selection. The boxes were 100 count s _1 wide 
in intensity and with a range of hardness ratio of 0.005. It 
was checked that the selections used were correct by over- 
laying the selected data on the full Z-track. In Fig. 2 it can 
be seen that the boxes are small in comparison with the Z- 
track as a whole so spectral fitting results relate accurately 
to a small part of the Z-track. It would clearly be inap- 
propriate to use large boxes extending along a significant 
fraction of the Z-track as appreciable spectral changes will 
take place within the box. However, although the boxes 
used in the present work are small, the count rates in the 
PCA and HEXTE are so large that the total count ac- 
cumulated in each box is very high as shown in Table 1 
allowing high quality spectral fitting. Accurate determina- 
tion of QPO frequencies was also possible, this requiring 
good count statistics as the fractional modulation is only 
a few percent. 

After initial tests, it was decided to use 11 spectra con- 
sisting of three spectra on the horizontal branch, one at the 
hard apex, four on the normal branch, one at the soft apex, 
and two in flaring, numbered as spectrum 1 (end of HB) 
to 11 (FB). An extra spectrum close to the peak of flaring 
(shown in Fig. 2) was found to have insufficient counts for 
good spectral fitting. Background spectra were generated 
for each selection. Both source and background data were 
deadtime corrected using a local facility pcadead. Pulse 
pileup correction has negligible effect on the spectra and 
so was not made. It was not necessary to regroup spectral 
channels to a specified minimum count to allow use of the 
X 2 statistic as the count in all channels was already high. 



A systematic error of 1% was added to each channel as 
usual for PCA data. 

HEXTE spectra and background files were produced 
for each selection using the same GTI files and deadtime- 
corrected using the ftool hxtlcurv. The auxiliary re- 
sponse file (arf) of May 2000, and the response matrix 
file (rmf) of March 1997 were used in spectral analysis. 
The rmf file was rebinned to match the actual number of 
HEXTE channels using RDDESCR and rbnrmf. For both 
the PCA and HEXTE, the source + background spectra 
were individually compared with background spectra and 
spectral fitting only carried out up to the energy at which 
these became equal, typically 30 kcV in the PCA and 40 
keV in HEXTE 

3.1. Spectral analysis 

The motivation of this work is to test the hypothesis that 
the extended ADC model (Sect. 1) can provide good fits to 
the spectra of the Z-track sources, and to test the physical 
model previously proposed to explain the Z-track source 
GX 340+0 (Church et al. 2006) in the case of GX5-1. 
Consequently, the extended ADC model was used in Xspec 
spectral fitting software in the form bb + CPL as discussed 
in Sect. 1. 

Preliminary fitting showed that the model provided 
good values of % 2 /d.o.f for all spectra. The residuals ex- 
hibited no evidence for an iron line; however, at the peak 
of flaring, there was evidence for the 10 keV emission fea- 
ture seen by Asai et al. (1994). This was seen clearly in 
the spectrum eventually discarded because of insufficient 
counts, but was also detectable in spectra 10 and 11 on the 
flaring branch at a significance of about 2a. When added 
to the fitting, the line had little effect on the results, for 
example, changing the neutron star blackbody radius by 
about 5%. Thus the model finally used in fitting all of the 
spectra omitted the line. 

In the Z-track sources, it is well-known that the 
Comptonized emission has a relatively low cut-off en- 
ergy of a few keV whereas lower luminosity sources have 
a higher cut-off energy which can approach 100 kcV. 
Broadband data such as from BeppoSAX or Suzaku allow 
measurement of the cut-off energy and the power law pho- 
ton index T. But in the Z-track sources including GX5-1, 
the low cut-off energy restricts the energy range available 
for determination of T and we have taken the approach 
previously adopted (e.g. Church et al. 2006) of fixing T 
at 1.7, a physically reasonable value (Shapiro et al. 1976). 
This procedure gave good quality spectral fits, and when 
a solution was obtained, it was found that T could be 
freed, but then remained close to 1.7. Extensive testing 
also revealed that the results obtained, i.e. the pattern of 
variation of parameters along the Z-track did not depend 
on the value of power law index. 

The lower energy limit of the PCA, usually set in spec- 
tral fitting at ~3 keV, does not allow very accurate deter- 
mination of the column density, and it is often the case 
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Fig. 3. Left: best-fit to the spectrum of the hard apex with residuals, right: the corresponding unfolded spectrum 



Table 2. Spectral fitting results: 90% confidence errors are shown. 



spectrum 




kT 
keV 


normsB 


km 


Eco 
keV 


normcPL 


X Vd.o.f. 


Horizontal Branch 

f 

2 

3 


6.1±0.5 
6.3±0.5 
6.5±0.4 


2.08±0.10 
1.98±0.10 
1.93±0.08 


7.6±0.7 
7.2±0.8 
7.4±0.8 


5.1±0.5 
5.4±0.6 
5.8±0.5 


5.9±0.3 
5.8±0.2 
5.6±0.2 


19.5±1.4 
21.2±1.5 
23.4±1.5 


32/58 
24/58 
46/58 


Hard Apex 
4 


6.8±0.5 


1.87±0.09 


7.2±0.1 


6.07±0.6 


5.6±0.2 


25.3±1.8 


38/58 


Normal Branch 

5 

6 

7 

8 


6.6±0.4 
6.5±0.5 
6.4±0.5 
6.2±0.5 


1.77±0.07 
1.70±0.10 
1.54±0.07 
1.41±0.06 


6.6±1.1 
5.6±0.6 
5.0±0.1 
4.8±0.2 


6.5±0.5 
6.5±0.8 
7.5±0.7 
8.9±0.8 


5.4±0.1 
5.1±0.1 
4.9±0.1 
4.7±0.1 


25.5±1.7 
26.9±0.8 
26.1±2.0 
24.4±2.3 


37/58 
41/58 
40/58 
38/58 


Soft Apex 
9 


5.7±0.7 


1.31±0.05 


6.0±0.2 


11 3+ 1 ' 4 


3.7±0.1 


25.1±3.3 


74/58 


Flaring Branch 

10 

11 


5.8±0.7 
5.9±0.7 


1.39±0.06 
1.46±0.05 


8.6±0.2 
11.6±0.1 


12. Oil. 
12.7±0.8 


3.6±0.2 
3.5±0.2 


25.6±3.8 
27.2±4.0 


48/58 
40/58 


Column densities are in units of 10 22 atom cm 2 ; the normalization of the blackbody is in units of 10 37 
for a distance of 10 kpc, the normalization of the cut-off power law is in units of photon cm~ 2 s -1 keV 


erg s _1 
- 1 at 1 



keV. 



that Nn has to be fixed at the Galactic value. In the 
present case, the Galactic value of ~ 0.94 x 10 22 atom 
cm~ 2 (Dickey & Lockman 1990) was tried, but the fits 
were unacceptable with, for example, x 2 /d.o.f of 288/59 
for spectrum 1. Christian & Swank (1997) also found a 
value higher than Galactic (2.54 x 10 22 atom cm -2 ). In the 
present case, with free JVh, the values varied between 5.7 
and 6.8 xlO 22 atom cm -2 along the Z-track and x 2 /d.o.f. 
became acceptable. Moreover, it was apparent that JVh 
changed systematically around the Z-track showing not 



only is the extra absorption intrinsic, but providing evi- 
dence for mass release or outflow within the system asso- 
ciated with the physical changes taking place. The same 
effect was seen in the case of RXTE observations of the 
Z-track source GX 340+0 (Church ct al. 2006). The best- 
fit to the spectrum of the hard apex is shown in Fig. 3 
as both the folded data (with residuals) and as unfolded 
data. Final fitting results are shown in Table 1 and in the 
following we discuss the behaviour of the blackbody and 
Comptonized emission. 
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Fig. 4. The neutron star blackbody emission: upper panel: the Fig. 5. Variation of emission component 1-30 keV luminosi- 



blackbody temperature; lower panel: the blackbody radius 



ties around the Z-track: upper panel: the blackbody luminosity; 
lower panel: the ADC Comptonized emission luminosity 



Firstly, we show in Fig. 4 the variation of the 
blackbody temperature /cTbb and blackbody radius Rbb 
around the Z-track as a function of the total 1-30 
keV luminosity Lxot for which a source distance of 9 kpc 
(Christian & Swank 1997) was assumed. The blackbody 
radius is defined by Lbb = 47ri?| B erT 4 where Lbb is 
the blackbody luminosity and a is Stefan's constant. The 
blackbody parameters change systematically along the Z, 
fcTsB having its lowest value of 1.3 keV at the soft apex 
where -Rbb is maximum at ll^J '\ km. Moving along the 
normal branch, kT-QB increases to 1.9 keV at the hard 
apex, and continues to increase on the horizontal branch, 
reaching 2.1 keV at the end of the Z-track. At the soft 
apex, the value of Rbb indicates that the whole neutron 
star is emitting, but decreases on the normal branch show- 
ing that the emission contracts presumably to an equato- 
rial belt. In Sect. 4.1 we discuss this measured decrease 
in emitting area on the neutron star showing that it is 
consistent with the effects of increased radiation pressure. 
Thus at the soft apex the source appears to be in an undis- 
turbed state with the neutron star emitting from its whole 
surface at its lowest temperature; we see below that the 
evidence indicates the mass accretion rate is minimum at 
this position. 

In Fig. 5 we show the variation of the luminosities of 
the individual emission components Lbb and Ladc : the 
neutron star blackbody and the Comptonized emission of 
the ADC. Firstly, it is clear that the Comptonized emis- 
sion is dominant, its luminosity being ten times larger than 
that of the blackbody. However, in flaring, the blackbody 
luminosity increases by a factor of two to ~40% of the to- 
tal at the peak of flaring. Moving away from the soft apex 
on the normal branch there is a strong increase of Ladc 



Table 3. Total luminosity, blackbody luminosity and ADC 
luminosity in the band 1-30 keV in units of erg s _1 . 



spectrum 


^Tot 




Lxdc 




(10 38 ) 


(10 37 ) 


(10 38 ) 


1 


6.30 


6.19±0.57 


5.68±0.44 


2 


6.74 


5.82±0.61 


6.16±0.44 


3 


7.21 


5.99±0.65 


6.62±0.41 


4 


7.67 


5.82±0.77 


7.08±0.50 


5 


7.51 


4.97±0.81 


6.98±0.46 


6 


7.44 


4.42±0.81 


6.98±0.60 


7 


7.06 


4.08±1.30 


6.66±0.54 


8 


6.62 


3.86±1.40 


6.24±0.61 


9 


5.52 


4.86±1.40 


5.04±0.68 


10 


5.74 


6.94±1.40 


5.06±0.75 


11 


6.07 


9.40±0.81 


5.12±0.75 



As is normal practice, we do not quote errors in the total 
luminosity. 



by 40% on the normal branch, with the PCA intensity in- 
creasing from ^10 000 to ^15 500 count s -1 . The increase 
in X-ray intensity is, of course, unquestionable. 

The significance of the ADC luminosity increase on 
the normal branch was assessed by fitting the variation 
of Ladc with Lxot (lower panel of Fig. 5) also shown in 
Table 3 using only the normal branch data. A linear fit 
gave a slope of 0.945±0.35, where the errors at 90% con- 
fidence strongly exclude a slope of zero. The probability 
that the slope is zero (i.e. testing the hypothesis that the 
luminosity is constant) is less than 0.006% so that this in- 
crease of ADC luminosity on the normal branch is highly 
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significant. Fig. 5 shows that the increase of X-ray inten- 
sity in climbing the normal branch is due to the increase 
in Ladc- We suggest that it is unlikely that Ladc could 
increase so substantially without an increase in mass ac- 
cretion rate M, and this is also indicated by the increase 
in neutron star blackbody temperature. In previous work, 
it was often assumed that the Z-track was produced by 
the mass accretion rate M varying monotonically in the 
direction HB - NB - FB (Sect. 1), although this conflicted 
with the known increase of X-ray intensity between the 
soft and hard apex. We propose here based on the mea- 
sured definite increase of Ladc that M increases on the 
normal branch between the soft apex and hard apex, in 
the opposite direction to that normally assumed. The be- 
haviour of the neutron star blackbody is more complicated 
consisting of a temperature increase but an area decrease 
such that the overall luminosity does not change substan- 
tially as discussed in Sect. 4.1 and so we should not expect 
Lbb to follow M in a simple way. 

The generally accepted idea that M increases mono- 
tonically around the Z-track explains flaring as a lumi- 
nosity increase due to high and presumably variable M. 
However, it is clear from Fig. 5 that Ladc is constant 
within errors on the FB indicating that M does not change 
in flaring. The blackbody luminosity however, clearly in- 
creases and we propose that flaring must consist of lu- 
minosity increase due to unstable thermonuclear burning 
on the surface of the neutron star, a suggestion that has 
not previously been made. During flaring, there is a rel- 
atively small increase of blackbody temperature, but the 
blackbody radius rises to ~13 km for spectrum 11, and 
would be even larger at the peak of flaring. Values of _Rbb 
substantially larger than the probable neutron star radius 
have been found previously in GX 5-1 (Christian & Swank 
1997), and a type of radius expansion known in X-ray 
bursts may be taking place, with the emitting region ex- 
panding beyond the normal surface of the neutron star. 

Figure 6 shows the variation of the ADC Comptonized 
emission parameters: the normalization and the cut-off en- 
ergy Eqo- Firstly, there is an increase of normalization 
moving on the NB away from the soft apex, while Eqo 
increases, these factors combining to give the observed in- 
crease of Ladc (Fig- 5). Secondly, comparison of Figs. 4 
and 6 shows that Eco is well correlated with /cXbb on 
the HB and NB, suggesting a link between the emission 
from the neutron star and the electron temperature of the 
ADC. For an ADC with high optical depth to electron 
scattering Eqo = 3fcT e (Petrucci et al. 2001) whereas for 
low optical depth Eqo — kT e . The present results allow 
a comparison of the neutron star blackbody temperature 
and Eco an d it is found that Eco/kT^B has an aver- 
age value of 2.7±0.5 implying that the plasma is optically 
thick and that the ADC electron temperature is close to 
the neutron star blackbody temperature. The relationship 
of fcTBB and kT e in LMXB was discussed by Balucihska- 
Church & Church (2005) who showed that in lower lu- 
minosity sources the ADC electron temperature kT e was 
substantially higher than kTsB so that there must be a 
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Fig. 6. The cut-off energy (upper panel) and normalization 
(lower panel) of the ADC Comptonized emission. 



presently-unknown heating process. However, in brighter 
sources (i ^ 2 x 10 37 erg s _1 ), these two temperatures 
were found to be consistently equal, providing evidence for 
thermal equilibrium between the neutron star and corona. 
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Fig. 7. Lower panel: The variation of lower frequency and up- 
per frequency kHz QPO detected on the upper normal and 
horizontal branches as a function of the total 1-30 keV lu- 
minosity (I/Tot) together with the difference Aw, top panel: 
variation of the frequency ratio V2/1/1 around the Z-track 
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3.2. Timing analysis 

We next examine the results of carrying out timing anal- 
ysis for the same selections of data along the Z-track 
as used in Sect. 3.1. Two types of data had been accu- 
mulated suitable for timing analysis up to kHz frequen- 
cies: event mode data of type E_125us_64M_24_ls and sin- 
gle bit data of types SB_125us_0_13_ls, SB.125us_14.17_ls 
and SB_125us_18_23_ls, i.e. one mode covering the energy 
band up to 20 keV stored in three files. Power density 
spectra were extracted with 125 /is time resolution from 
both event mode and single bit data into a single power 
spectrum for each selection along the Z-track. The energy 
range was restricted to the band 5-60 keV for event mode 
data with the i?XT#-specific ftool SEBITMASK specifying 
channels corresponding to these energies and applying this 
with feselect to the data files. Power spectra were then 
extracted using the Xronos facility powspec which essen- 
tially rcbinned all data into a single energy bin. A stan- 
dard Leahy normalization was applied. In order to carry 
out high quality fitting of the power spectra, we converted 
the spectra to a format that allowed fitting with the Xspec 
spectral fitting software, requiring changes to the header 
of the fits files and re-naming of the columns to corre- 
spond to frequency power and power error. Finally an 
output file was made suitable for use with Xspec using the 
ftool flx2xsp producing a file with pha extension as for 
an energy spectrum. This ftool also created a dummy in- 
strument response function as needed by Xspec consisting 
of a unitary matrix having no effect on the data. Deadtime 
correction was carried out and the spectra grouped to a 
minimum of 60 counts per bin as found suitable by trial 
and error. 

The power spectra were fitted within the range 150 
- 950 Hz, using a three-component model consisting of 
two Lorenzian lines for the kilohertz QPO plus a power 
law which was able to fit the continuum of the spectrum 
well within the restricted frequency range. As expected, 
kHz QPO were only detected on the horizontal branch 
of the Z-track and the upper part of the normal branch, 
that is, from a total of 5 of the 11 positions at which 
spectral fitting was carried out. Fig. 7 (lower panel) shows 
the best-fit frequencies obtained for the two QPO peaks 
in these spectra as a function of total luminosity (I/Tot) 
in the band 1- 30 keV obtained previously during spectral 
fitting. It can be seen that there is a monotonic systematic 
increase of the QPO frequency from the end of the HB to 
the point on the upper NB for both the upper and lower 
frequency QPO. Thus we find the same general trend of 
the kHz QPO frequencies as reported previously for GX 5- 
1 (Jonker et al. 2002) and in other Z-track sources (e.g. 
Jonkcr et al. 2000). The frequency difference v 2 - v\ is 
also shown in the lower panel and can be seen to have a 
minimum of ^280 Hz at the hard apex and a maximum of 
~380 Hz part way along the horizontal branch. In terms of 
the 90% confidence uncertainties in v 2 and v\ obtained in 
the fitting, typically 10 - 30 Hz, the variation is significant. 
The ratio of frequencies vijv\ also shown in the upper 



panel. Although the ratio passes through the value 1.5, 
i.e. the ratio 3:2 significant in the resonance model of QPO 
(Sect. 1), and interestingly this occurs at the hard apex, 
the ratio varies substantially around the Z-track from 1.4 
on the normal branch to 2.4 at the end of the horizontal 
branch, only equalling 1.5 in one part of the track. 
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Fig. 8. The orbital radii of the upper and lower frequency QPO 
as a function of luminosity 



We next assume that the upper frequency v 2 cor- 
responds to an orbital frequency in the inner disk and 
calculate the orbital radii for both kHz QPO assuming 
Newtonian dynamics in which the orbital frequency at a 
radius r in the disk is given by 



2-kv = 




where M is the mass of the neutron star assumed to be 
1.4 Mq. These radial positions are shown in Fig. 8. For 
the higher frequency QPO, the frequency v 2 for the power 
spectrum in the normal branch corresponds to a radial 
position of about 18 km, i.e. some distance away from the 
neutron star surface. 

As the source moves on to the horizontal branch, the 
frequency decrease corresponds to a radial position mov- 
ing outwards eventually reaching 25 km clearly suggesting 
that the source of oscillation moves further away from the 
neutron star. In terms of the spectral evolution along the 
Z-track, this is very suggestive, since both in the present 
work (Sect. 4.1; Fig. 9) and in our previous work on 
GX 340+0 (Church et al. 2006), the radiation pressure at 
the inner disk becomes very strong on the HB and is ex- 
pected to disrupt the inner disk removing the inner por- 
tions. Discussion of this will be postponed until after more 
detailed consideration of the spectral fitting results. 
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4. Discussion 

4.1. Spectral fitting results 

The results obtained in Sect. 3.1 show that the spectra 
of GX5-1 are well- fitted by the extended ADC model, 
and that the evolution of spectral parameters around the 
Z-track suggest a relatively simple explanation of the Z- 
track as discussed below. The increase of ADC luminosity 
between the soft apex and hard apex provides evidence 
that M increases in this direction. The constancy of the 
ADC luminosity on the flaring branch combined with an 
increasing Lbb shows that flaring must consist of unstable 
nuclear burning on the neutron star. Both of these conclu- 
sions disagree with the often-held view that M increases 
monotonically in the direction HB - NB - FB. Moreover, 
the spectral evolution of GX5-1 around the Z-track is 
essentially the same as in the similar source GX 340+0 
(Church et al. 2006) so that the same explanation works 
for both sources. We discuss below effects taking place 
on the normal and horizontal branches and on the flaring 
branch in more detail. 



4.2. The effects of radiation pressure 

Moving away from the soft apex where the neutron star 
blackbody temperature is lowest and the blackbody radius 
indicates that the whole star is emitting, we find a sub- 
stantial increase in fcTee , so that there is a large increase 
in the radiation pressure ~ T 4 . 
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Fig. 9. The ratio of neutron star emitted flux to the Eddington 
flux as a function of position on the Z-track. 



At the same time, the blackbody radius Rbb decreases 
to a value of ~5 km at the HB end of the Z-track so 



that the emission is from a reduced area on the neutron 
star, presumably an equatorial belt, ft would be normal 
to describe the source in terms of its luminosity compared 
with the Eddington luminosity which is, of course, derived 
on the assumption of radial and not disk accretion and 
assuming accretion over 47T steradians. In the present case, 
where only parts of the neutron star emit, it is instructive 
to compare the flux / emitted by unit emitting area on the 
neutron star with the Eddington flux LEdd/47ri? 2 , where 
R is the neutron star radius assumed to be 10 km, so that 
/Edd = 1-4 x 10 25 erg cm~ 2 s _1 . When only parts of the 
neutron star are emitting, the flux / will exceed the mean 
X-ray flux deduced from the measured luminosity and the 
source distance. Figure 9 shows the ratio of / to /Edd- The 
ratio has its lowest value at the soft apex of ~0.2 showing 
that radiation pressure is weak, rising to unity at the hard 
apex and continuing to rise to the end of the Z-track. 
Thus the radiation pressure of the emitting band becomes 
strong at the hard apex and horizontal branch such that 
the effects on the accretion flow will be strong. At the 
Eddington flux, the radiation pressure force balances the 
inwards gravitational force, both flux and gravitational 
force decreasing as 1/r 2 as in the standard derivation of 
the Eddington limit, except that this does not apply over 
the whole neutron star surface and the effects of radiation 
pressure will depend on which direction is considered. 

Especially in luminous LMXB, the inner accretion 
disk is radiatively supported by the radiation pressure of 
the disk itself having an equilibrium half-height iJ cq = 
Sa^M /8irm p c, where ctt is the Thomson cross section, 
m p is the proton mass and c the velocity of light (Frank 
et al. 2002). The height of the disk rises rapidly from zero 
to H eq in a radial distance ~10 km from the surface of 
the neutron star, and for the hard apex in the present 
data, H cq is 65 km where M is derived from the total 1 
- 30 keV luminosity. Thus the inner disk rises high above 
the neutron star, and we propose that the strong radia- 
tion pressure at the hard apex and on the HB severely 
disrupts the inner disk blowing large parts of it away. At 
any radial and vertical position in the unperturbed inner 
disk (r, z), the vertical force due to radiation pressure 
plus gas pressure is balanced by the vertically downwards 
component of the gravitational force. Radiation pressure 
from the neutron star acts on a line from the star to this 
point (r, z) offsetting the gravitational force. If this off- 
set reduces the effective gravitational force close to zero, 
the equilibrium height of the disk would tend to infinity. 
Moreover, if we assume that M is maximum on this part 
of the Z-track, the radiation pressure of the disk itself will 
also be increased moving the disk surface upwards. 

Clearly the combined effects of the radiation pressure 
of the neutron star and that of the inner accretion disk 
will be complicated and modelling should be carried out. 
That is beyond the scope of the present work, but we show 
in Fig. 10 an approximate schematic view of the expected 
effect. The unperturbed profile of the disk is shown as a 
dashed line, and the actual profile of the inner disk by 
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Fig. 10. Schematic view of the inner, radiatively supported 
inner accretion disk disrupted by radiation pressure; the un- 
perturbed disk is shown as a dashed line. The radial profile of 
height H(r) meets the unperturbed profile at radial position 
r un and it is proposed in the present work that the higher fre- 
quency kHz QPO consists of an oscillation taking place in the 
neighbourhood of r un 



a full line. The disk height H(r) will revert to the un- 
perturbed form of the radiatively-supportcd disk at some 
radial position r un , and this position represents an effec- 
tive inner edge of the disk, since at smaller radial positions 
there may be a residual disk, but of substantially reduced 
height. 

The increase of column density found in spectral fit- 
ting on the normal branch supports the release of mass 
within the system on this part of the Z-track. Our mea- 
sured increase of absorption indicates the absorber is spa- 
tially extended so that both emission components will be 
affected. It is clear that most of this material would not 
be fully ionized since the ionization parameter £ = L/nr 2 
decreases rapidly with distance and a £ of 10 6 is required 
for there to be complete ionization (e.g. Makishima 1986). 
Moreover there would be a spatial variation of £ so that 
material located at vertical positions not directly facing 
the emitting part of the neutron star would be less ion- 
ized. 

The localized radiation pressure of the neutron star as 
measured by ///Edd rises by almost an order of magnitude 
from the soft apex to the end of the HB. Radio emission is 
non-existent on the FB and soft apex and strongest at the 
hard apex and HB. Thus in GX5-1 as found previously 
in GX 340+0 there is a very good correlation between ob- 
servation of radio and / / /Edd and we propose that strong 
radiation pressure is a necessary, but not necessarily suf- 
ficient condition for jet formation. 

We also comment on the observed decrease of black- 
body radius, i.e. emitting area on the neutron star in as- 
cending the normal branch, noting that this is consistent 
with previous work. The blackbody radius obtained from 
spectral fitting Rbb is related to the half-height of the 
emitting equatorial strip on the neutron star h by 4 ir i? BB 
= Air i?NS h, where i?NS is the radius of the neutron star, 
assumed to be 10 km. In a survey of LMXB with Asca and 



BeppoSAX, Church & Balucinska-Church (2001) made a 
survey comparing the blackbody luminosity (Lbb) in the 
survey sources with the Comptonized luminosity (£adc) 
in order to investigate the well-known discrepancy that 
Lbb in most LMXB is many times less than the fraction 
of the total luminosity expected (^50%) on the basis of 
simple theory. It was found that the half-height of the 
equatorial emitter on the neutron star h was in the lower 
luminosity LMXB substantially less than the neutron star 
radius. However, in all sources h followed the simple re- 
lation h = H, where H is the half- height of the inner 
disk, valid over 3 decades of source luminosity. A proba- 
ble explanation of this geometric equality is the model of 
accretion flow spreading on the neutron star (Inogamov 
& Sunyaev 1999), which was shown to agree reasonably 
with the survey results (Church et al. 2002). Thus in lower 
luminosity sources, the spreading does not extend so far 
vertically, and the emitting region is smaller. 

We have suggested that in the Z-track sources parts of 
the inner accretion disk are blown away by high radiation 
pressure presumably leaving a residual innermost disk of 
reduced height. Thus a reduced value of disk height H is 
expected to lead to a reduced emitter size on the neutron 
star, and conversely, the observed decrease of i?BB and h 
is evidence that the height of the inner disk H is reduced 
from more than 50 km to ~3 km indicating disruption of 
the inner disk. 

4.3. Flaring 

The constancy of Ladc on the flaring branch and the in- 
crease of L B b lead to the conclusion that flaring must be 
unstable nuclear burning on the neutron star, since if M 
increased as usually assumed we would definitely expect 
Ladc to increase substantially as it does on the normal 
branch. Further support for this conclusion comes from 
comparing with the theory of unstable nuclear burning. 
Unstable burning is expected at various mass accretion 
rates (Fujimoto et al. 1981; Fushiki & Lamb 1987; Bildsten 
1998; Schatz et al. 1999), the various regimes of stable 
and unstable burning depending on the value of to, the 
mass accretion rate per unit area of the neutron star. For 
5.0 x 10 3 < to < 1.3 x 10 5 g cm~ 2 s _1 there is unstable He 
burning in a mixed H/He environment, while for to > tost 
= 1.3 x 10 5 g cm~ 2 s _1 , burning is always stable (Bildsten 
1998), this value having an estimated uncertainty of 30% 
(giving a range of about 0.9 - 1.7xl0 5 g cm" 2 s _1 ). In 
deriving values of to, we use actual values of the emitting 
area given by the blackbody radius obtained in spectral 
fitting, i.e. to = M /(4-7ri?| B ). At the soft apex, we find to 
= 1.8 ± 0.2 x 10 5 g cm~ 2 s _1 consistent with the critical 
value tost, and so it is expected theoretically that burn- 
ing is stable on the normal branch but unstable on the 
flaring branch, the soft apex being the point at which the 
transition between these occurs. 

If we progress along the FB away from the soft apex, 
the value of to derived from spectral fitting decreases, 
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so that once started, unstable burning remains unstable. 
Moving in the opposite direction away from the soft apex, 
however, m rises substantially, to 3.6 ± 0.5 x 10 5 g cm~ 2 
s^ 1 for the first spectrum on the NB, increasing eventu- 
ally to 10 x 10 5 g cur 2 s" 1 at the end of the HB. This 
is due to the decrease of blackbody emitting area and the 
increase of M and shows that conditions on the neutron 
star are well into the regime of stable burning. 

It would be desirable to find other observational tests 
of whether flaring consists of nuclear burning and it would 
be interesting to derive an a-parameter for flaring in GX 5- 
1, since the agreement of measured a values with the- 
ory was the essential proof that bursting was nuclear in 
nature, a is defined as the ratio of accretion luminosity 
to burst luminosity evaluated in bursting by integrating 
the non-burst luminosity between bursts, and the lumi- 
nosity within a single burst. For a single burst, the time 
needed for the accretion flow to accumulate and to settle 
on the neutron star increasing in density towards igni- 
tion conditions is well-defined as the time between bursts. 
Theoretical values are obtained by dividing the accretion 
luminosity per nuclcon (GM/Rns by the nuclear energy 
release per nucleon from which a ^30 for H-burning and 
more than 100 for He burning. However, a detailed inves- 
tigation of flaring in the Z-track sources in the context of 
nuclear burning has not been made and the physics of the 
accumulation of mass and settling on the star has not been 
investigated and it is unclear what timescale is involved in 
the matter involved coming to ignition conditions. Thus 
consideration of the a-parameter is beyond the scope of 
the present work. 

4.4. The kHz QPO 

We are now able to discuss the QPO results in terms of the 
spectral fitting results, and we firstly take the QPO vari- 
ation along the Z-track shown in Fig. 8 as the Keplerian 
radius of the oscillation versus Lxot and re-draw this as 
Keplerian radius as a function of ///Edd as shown in Fig. 
11. 

A striking feature can be seen that the major change in 
Keplerian radius and thus in QPO frequency takes place 
when ///Edd equals unity. This result is unlikely to be co- 
incidental and suggests that the strong radiation pressure 
causes the variation in QPO frequency along the Z-track. 
Moving from the hard apex to the end of the Z-track, 
///Edd increases from 1.0 to 1.5 and the Keplerian ra- 
dius of the higher frequency QPO increases as the fre- 
quency decreases. The most obvious explanation is that 
this radiation pressure is indeed disrupting the inner disk 
so that the oscillations comprising the higher frequency 
QPO move outwards but staying on the inner edge or 
close to it, and this can explain the change of frequency. 
It is clear that the lower frequency QPO responds to the 
change in in some way but the results do not provide 
a clear indication of the mechanism of this. However, the 
result shown in Fig. 11 does suggest an alternative to ex- 



isting ideas regarding possible truncation of the inner disk. 
Much theoretical work suggests that the accretion disk is 
in contact with the neutron star, e.g. the Inogamov & 
Sunyaev model (1999) of accretion flow spreading on the 
surface of the neutron star, while other work considers the 
possibility that the disk is truncated at a specific radius, 
either at an innermost stable orbit, or at an Alfven radius 
assumed external to the neutron star by capture onto the 
weak magnetic field of the neutron star or at a sonic point. 
In the sonic point model (Miller ct al. 1998), part of the 
accretion flow can remain in the disk until it impacts on 
the neutron star. The present work suggests that the disk 
may be truncated, but only in part of the Z-track, and 
that the mechanism may be radiation pressure. 

Rossi-XTE PCA GX 5-1 
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Fig. 11. Variation of Keplerian radius of the kHz QPOs as 
a function of radiation pressure of the neutron star, i.e. the 
emitted flux of the surface in units of the Eddington flux / E dd 
demonstrating a striking correlation between the QPO radial 
position of origin and high radiation pressure 



5. Conclusions 

We have shown that spectra selected along the Z-track 
in GX5-1 are well fitted using the extended ADC model, 
and that the spectral evolution is very similar to that we 
found in GX 340+0 (Church et al. 2006). Thus the phys- 
ical model proposed to explain the spectral evolution in 
that source is also able to explain GX5-1. In this model, 
the soft apex has the lowest value of M and the strong 
physical changes on the normal branch are primarily due 
to an increase of M and of neutron star temperature. The 
increase of temperature means that the radiation pressure 
must be high in the neighbourhood of the neutron star. 
We propose that this super-Eddington radiation pressure 
combined with inner disk radiation pressure disrupts the 
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inner disk diverting accretion flow into the vertical direc- 
tion leading to jet formation. Thus high radiation pressure 
is a necessary but not necessarily sufficient condition for 
jet formation. On the horizontal branch, the decrease of 
intensity and Lxdc indicates that M decreases; however 
the blackbody temperature continues to increase and there 
are various ways of explaining this, such as M decreasing 
in the disk while an enhanced accretion flow is still af- 
fecting the neutron star. Evidence is presented that the 
flaring branch consists of unstable nuclear burning on the 
neutron star given the constant Labc and the increase of 
Lbb, and the good agreement between m at the soft apex 
and the theoretical ttist below which burning is unstable. 
However, we should add the caution that while we now 
have evidence that the flaring branch can be explained 
in this way in the three Cyg-like Z-track sources, i.e. in 
GX 340+0 (Church ct al. 2006), in GX5-1 (present pa- 
per) and in CygX-2 (Church et al. 2007), we do not claim 
that this is so in the other group of Z-track sources: the 
Sco-like sources. Preliminary evidence suggests that these 
are similar on the NB and HB, but that the FB may be 
different. 

Radio emission is detected at the hard apex and hor- 
izontal branch showing that jets are present (Penninx 
1989; Hjellming & Han 1995; Berendsen et al. 2000) where 
we find that radiation pressure is strongest and is non- 
existent on the FB and soft apex. It was, of course, pre- 
viously suggested that radiation pressure may be impor- 
tant in jet formation (Bisnovatyi-Kogan & Blinnikov 1977; 
Begelman & Rees 1984) and that the conical opening in 
the inner accretion disk may perform some degree of col- 
limation (Lynden-Bell 1978). 

The present results do not support the standard view 
that M increases monotonically in the direction HB - NB 
- FB, but have M constant in flaring and increasing from 
the soft apex to hard apex. Moreover, in the standard 
view the source becomes super-Eddington at the soft apex 
whereas the present work indicates the opposite, that this 
happens at the hard apex. 

By combining timing and spectral analysis we have 
shown that the variation of the kHz QPO along the Z- 
track is highly correlated with ///Edd, and the major 
change in the QPO frequencies takes place when this ra- 
tio equals unity, and the radiation pressure becomes large. 
This suggests that the higher frequency kHz QPO is an os- 
cillation at the inner edge of the disk, i.e. the edge formed 
by the action of radiation pressure removing a large sec- 
tion of the unperturbed disk, and that the edge moves 
to larger radial positions when the radiation pressure in- 
creases. A major feature of kHz QPO is that they are 
essentially a feature of the horizontal branch, fading away 
at the upper end of the normal branch and this is consis- 
tent with our model for the Z-track as the high radiation 
pressure of the neutron star revealed by spectral fitting 
is also a feature of the horizontal branch. Fig. 11 shows 
that when ///Edd has fallen to ^0.5, kHz QPO are no 
longer seen. Thus there is an indication that the kHz QPO 
only exist when part of the inner disk is blown away. The 



smallest Keplerian radius we measure has a value of 18 
km on the upper normal branch and if QPO continued to 
exist on the normal branch, the Keplerian radius would 
decrease further. However, as the radiation pressure falls 
there would be no truncation of the disk at an edge acting 
as a site for the oscillation, and this is consistent with the 
lack of kHz QPO on the lower normal branch and flaring 
branch. 

Thus there is a simple explanation for the known oc- 
currence of the higher frequency kHz QPO at a preferred 
radial position in the disk: that this is the disk edge, 
and the results suggest that the variation of frequency 
is caused by change in position of the edge. 

The physics of the accretion flow in the inner disk is 
very uncertain (e.g. van der Klis 2000), and probably the 
most important question is whether the disk in general 
in LMXB meets the neutron star or is truncated at an 
innermost stable orbit, an Alfvcn radius or a sonic point, 
and this is relevant to the QPO models which invoke trun- 
cation. However, it is difficult to say whether truncation 
of the disk due to relativistic, magnetospheric or other 
effects actually takes place. The present work presents a 
simple observational result on the higher frequency QPO 
and does not directly say anything about the lower fre- 
quency oscillation, and it is not appropriate to comment 
on the detailed physics of the QPO models here. However, 
the present work indicates that the inner disk is not in 
general truncated, but can become truncated on the up- 
per NB and HB. Disk truncation is due to high radiation 
pressure and we do not need to invoke truncation because 
of an innermost stable orbit, or an Alfven radius external 
to the neutron star or a sonic point. 
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